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Abstract 

We establish the validity of the heavy traffic steady-state approximation for a single 
server queue, operating under the FIFO service discipline, in which each customer 
abandons the system if his waiting time exceeds his generally-distributed patience 
time. This follows from early results of Kingman when the loading factor approaches 
one from below, but has not been shown in more generality. We prove the convergence 
of the steady-state distributions of the offered waiting time process and their moments 
both under the assumption that the hazard rate of the abandonment distribution is 
scaled and that it is not scaled. As a consequence, we establish the limit behavior of 
the steady-state abandonment probability and mean queue-length. 
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1 Introduction 

There have been several papers that have studied queueing systems with customer aban¬ 
donments. The knowledge of the long time asymptotic behavior (stability) of such systems 
is of great importance in practice; however, excepting special cases, the models of inter¬ 
est are too complex to analyze directly. The heavy traffic limits often provide tractable, 
parsimonious, and yet very meaningful approximations. In this paper, we consider a sin¬ 
gle server queue, operating under the FIFO service discipline, with generally-distributed 
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patience time (the GI/GI/I + GI queue). To the best of our knowledge, the only re¬ 
sults available in the literature in the GI / GI /I + GI setting are process-level convergence 
results. Our aim is to fill this gap by establishing the convergence of the steady-state distri¬ 
butions/moments of the offered waiting time process for the GI/GI/1 + GI queue in heavy 
traffic. As a consequence, we derive the limit behavior of the steady-state abandonment 
probability and also the mean steady-state queue-length. These results are necessary to 
rigorously solve steady-state performance optimization problems in the heavy traffic limit. 


Our asymptotic analysis relies heavily on past work that has developed heavy traf¬ 
fic approximations for the GI/GI/I + GI queue , using the of f ered w aiting time process. 
The offered waiting time process, introduced in iBacelli et al.l (jl984l i. tracks the amount 
of time an infinitely patient customer must wait for service. Its heavy traffic limit when 
the abandonment dis t ributi on is left unsealed is a reflected Ornstein-Uhlenbeck process 
(see Ward and Glvnn ( 20051 )). and its heavy traffic limit when the aband onment distribu¬ 
tion is scaled through its hazard rate is a reflected nonlinear diffusion (see iReed and Ward 
( 20081 ) i. However, those results are not enough to conclude that the steady-state distribu¬ 
tion of the offered waiting time process converges, which is the key to establishing the limit 
behavior of the steady- state abandonme i it pro bability and mean queue-length. Those lim¬ 
its were conjectured in Reed and Ward ( 2008l l. and shown through simulation to provide 
good approximations. However, the proof of those limits was left as an open question. 

When the system loading factor is less than one, since th e GI/GI/ I + GI queu e is 


dominated by the GI/GI/I queue, the much earlier results of iKingmanI (|l96ll . Il962l i for 


the GI/GI/I queue can be used to establish the heavy traffic steady-state approximation 
for the GI/GI/I + GI queue. The difficulty arises because, in contrast to the GI/GI/I 
queue, the GI/GI/I + GI qu eue has a steady-sta te distribution when the system loading 
factor equals or exceeds 1 (see Bacelli et al. ( 1984I B. Furthermore, the results in the afore¬ 
mentioned papers only provide the convergence of the steady-state distribution, but not 
its moments. 

An informed reader would recall Gamarnik and Zeevi ( 2006l i and Budhiraia and Lee 
( 20091 ). which establish the validity of the heavy traffic steady-state approximation for a 
gener alized Jackson network, wit hout customer abandonment. The proof of the former 
paper Gamarnik and Zeevi ( 200(tI ) relied on certain exponential integrability assumptions 
on the primitives of the network and as a result a f orm of exponential ergodicity was estab¬ 
lished. The latter paper Budhiraja and Lee ( 20091 ) provided an alternative proof assnming 
the weaker square integrability conditions that are commonly used in he avy traffic analysis. 
Our a nalysis is inspired by the methodology developed in the latter work iBudhiraia and Lee 
( 2 OO 9 I ) and we expand their approach to accommodate various customer abandonment set¬ 
tings. The key differences are that (a) the presence of customer abandonments means we 
must also use the properties of a (nonlinear) generalized regulator map in our analysis; 
(b) we have to worry about the difference between the offered waiting time process, which 
does not include customers that will eventually abandon, and the queue-length process, 
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that does; (c) the scaling of the hazard rate distribution leads to a connection between 
the tail behavior of the limiting hazard function and the number of moments that must be 
assumed on the inter-arrival and service times to ensure the convergence of the steady-state 
moments (which is not required to obtain distributional convergence only). 


In comparison to results for many-server queues, the process-level convergence re¬ 
sult for t he GI/GI/N + GI queue in t h e qua lity-and-efficiency-driven regime was estab- 
lis hed in iMandelbaum and Momcilovid (120121 ) when the hazard rate is not scaled, and 


m 


Reed and Tezcan ( 20121 ). under the assumption of exponential service times, when the 


hazard rate is scaled. Neither paper establishes the convergence of the steady-state distri¬ 
butions. That convergence is shown under the assumption that the ab andonmen t distr i- 
bution is exponential and the service time distribution is phase type in Dai et al. ( 20141 ). 
The question remains open for the fully general GI/GI/N + GI setting. 


The remainder of this paper is organized as follows. In Section [21 we set up the model 
assumptions and recall the known process-level convergence results for the GI/GI/! + GI 
queue both when the hazard rate distribution is not scaled and when it is. In Section |3l 
we state our main result, that gives the convergence of the steady-state distribution of the 
offered waiting time process, and its moments. As a corollary to this result, we obtain the 
convergence of the steady-state abandonment probability and mean queue-length. Section 
[4] shows how to obtain bounds on the moments of the scaled state process that are uniform 
in the heavy traffic scaling parameter (n). We use those bounds in Section [5] to prove our 
main result. 


2 The Model and Known Results 


The GI/GI/1+GI model having FIFO service is built from three independent i.i.d. se¬ 
quences of non-negative random variables {ui,i > 1}, {u*, i > I}, and {di,i > I}, that are 
defined on a common probability space (D,J^, iP), and for which IE[ui] = iEfui] = 1 and 
IE[ul] < oo,IE[vf] < oo. 

We consider a sequence of systems indexed by n > I in which the arrival rates become 
large and service times small. (By a convention, we will superscript any process or quantity 
associated with the n-th system by n.) More specifically, assume a sequence of arrival rates 
indexed by n > 1 defined as 

(Al) A” = nX for some A > 0. 


The z-th arrival to the system having arrival rate A"' occurs at time 


,n _ //2_ 

b \n’ 


1=1 
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and has service time 


where /i” is a service rate in the n-th system. A customer abandons without receiving 
service if processing does not begin by time + We let F be the cumulative distribution 
function of di, and assume it is proper (that is, F{x) —>■ 1 as x —>■ oo). 

We impose the standard heavy traffic assumption. 

(A2) y/n —)• d e (— 00 , 00 ) as n ^ 00 . 

For example, when 


fF = nX — y/nO + o{y/n) for some 0 € (— 00 , 00 ), 

so that the mean service times become short as the demand becomes large, (A2) holds. 
The arrival and service rates are approximately equal for large n because (Al) and (A2) 
imply 

iF A” 

- > X as re —)• 00 and — = A for every re. 

re re 

Finally, the following assumption will be required to ensure each system in the sequence 
is stable. 


(A3) The inter-arrival and service time distributions are such that 

sup{x ; lP{vi < /r”x) = 0} — inf{x : fP(rei < A”x) = 1} < 0 for each re > 1. 

The Offered Waiting Time and Qneue-Length Processes 

The offered waiting time process, first given in Bacelli et al. ( 1984l i. tracks the amount 
of time an incoming customer at time t has to wait for service. That time depends only 
upon the service times of the non-abandoning customers already waiting in the queue, that 
is, those waiting customers whose abandonment time upon arrival exceeds their waiting 
time. For t > 0, the offered waiting time process having initial state F”(0) = Xq has the 
evolution equation 


A^{t) 

V^it) = a;o ^ - 

i=i 

where 

A”'(t) = max{i : < t}. (2) 




( 1 ) 
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The quantity V^{t) can also be interpreted as the time needed to empty the system from 
time t onwards if there are no arrivals after time t, and hence it is also known as the 
workload at time t. 


The queue-length process Q^{t) represents the number of customers that are in the 
system at time t > 0, either waiting or with the server. In contrast to V”, includes 
customers that will eventually abandon but have not yet done so. Fortunately, the number 
of such customers is small enough that they can be safely ignored in our asymptotic regime, 
as will be seen. The implication for our analysis is that we can focus on V"", and obtain 
results on Q"’ from results on V"'. 


Th e steady-state distributi ons of V"" and exist uniquely from Bacelli et al. ( 19841 ) 


and/or Lillo and Martini ( 200lh . We let (oo) and Q^{oo) be the random variables having 
the respective steady-state distributions. 


Conventional knowledge states that the queue-size is of order of the square-root of the 
arrival rate, and so a sample path version of Little’s law known as the snapshot principle 
suggests that 




Q^jt) 

A*" 


1 1 

, for t > 0. 

y/n X 


The implication of the offered waiting time process becoming small is that it is only the 
most impatient customers that abandon. More specifically, assume F'(0) is well-defined 
and the initial condition satisfies y/nx^ —)• a: as n —>■ oo. Let {V,L) be the unique solution 
to the reflected stochastic differential equation 


V{t) =x + aW{t) + jt- F'(0) /o V(s)ds + L(t) > 0 

subject to: L is non-decreasing, has L(0) = 0 and V(s)dL(s) = 0, 


where {W{t) ■ t > 0} denotes the one-dimensional standard Brownian motion, and the 
infinitesimal variance parameter is 


cr^ = A var(Mi) -|- A ^ var(ui) 


(See t he paragraph below (I19p for a derivation.) Then, Theorems 1 and 3 in IWard and Glvnn 
( 20031 ) show that 


y/nV^ ^ V and —= ^ AlA, in D{]R+,M) as n —>■ oo, 
y/n 


(4) 


where D{]Rj^,]R) is the space of right-continuo us function s with left limits, endowed with 
the Skorokhod Ji-topology (see, for example, Billingsley (Ii 999I B. and the symbol 
stands for the weak convergence. 

The convergence in (j3|) motivates approximating the scaled steady-state distributions 
for C"' and Q”, and their steady-state moments, using the steady-state distribution of V, 
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a nd its moments. Thi s is co nvenient because there is the analytic expression for ^^[^(oo)] 
Browne and Whitt ( 1995l h which gives 


in 


]E[V{oo)] =x + 


e a -e I 2 


F'(0) 


where hz{-) is the hazard rate function of the standard normal distribution. When the 
abandonment distribution is exponential with rate 7 , F'{Q) = 7 , and the parameter 7 fully 
characterizes the abandonment distribution. Otherwise, there is the potentially undesirable 
feature that the approximation depends only on the value of the abandonment density at 

0 . 

Incorporating the Entire Abandonment Distribntion 

We apply the hazard rate scaling in Reed and Ward ( 2008l h which captures the full 
abandonment distribution. To do this, we must allow the abandonment distribution to 
change with n. Specifically, we replace the sequence {di^i > 1} with {d”,i > 1}, and let 
be the cumulative distribution of d”, and /i” its associated hazard rate function. Note 
that when E” has support on [ 0 , 00 ), then the relationship 


F^{x) = 1 — exp J hF{u)dv^ , a; > 0 , 


(5) 


holds, and we assume that is a nonnegative and continuous function on [ 0 , 00 ) 


To have an approximation in which the entire abandonment distribution appears, it 
is necessary to modify the diffusion in ([3]). For that, we begin by considering a customer 
arrival at time s > 0, that hnds the offered waiting time to be R”(s). The probability that 
would-be customer abandons is 


F"(R”(s)) = 1 - exp - 




hP'{u)du 


from (j5|). The recollection that waiting times are of order \l\pn suggests the change of 
variable u = ^/nv, because then 

_1 fVnV'^is) ^ f y 


/ _1 fy/nV (s) / V \ 

F^{V^{s)) = 1 - exp I — j dv 1 , 


and we expect the upper limit of integration to be convergent. Next, a Taylor series 
expansion suggests the approximation 




^ rV^V’^is) 


y/n 


L 


h^[ — \ dv. 

n. 
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The assumption that 


/i”(x) = h{-\/nx) 


( 6 ) 


for some given hazard rate function h implies that the scaled instantaneous loss rate for 
these would-be customers is 


(F”(s)) 



(7) 


where F(s) is the assumed limit of ^/nV'^{s). That scaled instantaneous loss rate appears 
as the instantaneous drift term F'{Qi)V{s) in (l3|). Replacing that term with the right-hand 
side of d?]) leads to the revised stochastic differential equation 


V{t) = X + aW(t) + jt — Jq h(u)du ds + L{t) > 0 

subject to: L is non-decreasing, has L(0) = 0 and V(s)dL(s) = 0. 


( 8 ) 


We note that the existence of a unique stationary distribution l/(oo) is always guaranteed. 
This is because h(u) > 0 is a hazard rate function, it satisfies h{u)du = oo, and 
therefore, trivially there exists a constant zq > 0 snch that Jq h{u)du > 9 for all z > zq. 


Theorems 5.1 and 6.1 in iReed and WardI ((20081) provide rigourous snpport (i.e., weak 
convergence as in @) for the approximation of y/nV^ by R in ([8]). This requires the 
following additional assumption. 


(A4) The abandonment distribution function F” and its associated hazard rate function 
h” satisfy 

F"(R"(s)) = 1-exp ( ^ y ^ j 1 and h^{x) = h{y/7lx), for all x > 0. 

Remark 1. (Connecting the Diffusion Approximations (0) and (0)^ The diffusion defined 
by (0) is in some se nse a refinement of the one defined by To see this (following 

intuition provided in Reed and TezcaA i '20lA) for the GI/M/N + GI model), perform a 
Taylor series expansion of h{u) about 0 as follows 


OO i 

h(u) = /.(O) + 

Then, since h{0) = F'(0), it follows that the linear portion of the infinitesimal drift in 
m arises from the lowest order term of the above Taylor series expansion. Hence the 
simpler diffusion approximation can only he expected to perform well when /i(0) is large 
compared to h^'^\0). This is true for an exponential distribution (since the hazard rate is 
constant). This is not true for the many Gamma distributions whose density equals zero 
at the origin. In such cases, the diffusion approximation in ^ is better able to capture 
the effects of customer abandonments by taking into account the entire behavior of the 
abandonment distribution. 
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t —)• oo 


y/nV^{t). 


,y/nV^{oo) 


n —>■ oo 


V{t) 


t —)• oo 


n —)> oo 


V (oo) 


Figure 1: A graphical representation of the limit interchange. 


3 The Steady-State Convergence Result 

We establish a limit interchange result, represented graphically in Figure [TJ Specihcally, 
for V dehned by (l3|) or V dehned by ([8]), it is known that 

V =^/nV'^ ^ V in D{lRj^,lR) as n —oo. 

However, the convergence of the associated steady-state distributions and steady-state 
moments is not known. For the convergence of the steady-state moments, we require the 
following slightly higher moment assumption on the primitives. 

(A5) Assume that 

E[u\ -|- < oo for some q > 2. 

Theorem 1. Let V be defined by ^ z/(Al)-(A3) hold and by i/(Al)--(A4) hold. 

(a) As n —^ oo, H"'(oo) ^ H(oo). 

(b) If also (A5) holds, then also 

iF[(H"(oo)n ^ iF[CF(oo)n, asn^oo (9) 

holds for any m G {0,q — 1), where q is as in (A5). 

The proof of Theorem [1] is found in Section O after appropriate technical machinery is 
developed in Section [H 

As a consequence of Theorem [H we are able to establish the limit behavior of the scaled 
steady-state abandonment probability and of the scaled queue-length. For the hazard rate 
scaling case (A4), we impose the following conditions on the hazard rate and moments of 
the model primitives. 
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(A5') Assume that 


h{u) < K{1 + v!") for some I > 0 and K > 0, for all tt > 0, 


and that 

lE[ul + < oo for some q G {2 + l,oo). 

Corollary 1. Let l/(oo) be a random variable possessing the steady-state distribution of 
V in Then, under (A1)~(A3) and (A5), we have 


V^P:^F'iO)]E[Vioo)] 


and 



A-E[l/(oo)], as n —>■ oo. (10) 


Let F(oo) be a random variable possessing the steady-state distribution ofV in (E\)- Then, 
under (A1)-(A4) and (A5'), we have 





and 



X]E[V{oo)], as n —)■ oo. 


( 11 ) 


Proof of Corollary We consider the case that (A1)-(A5) hold, and V is defined by ([8|). 
The case when (A4) is not assumed is a special case of the former by setting the abandon¬ 
ment distributions -F” = F and h{u) = F'(0) in what follows. 

To prove the first assertion in (|lll) . notice that 


gn{x) = ^/n (l - exp { 


X, 


as n —)• oo, uniformly on compact sets. This in turn implies that gn(xn) x whenever 
Xn —^ X. From Theor e m [Ha) and (gener alized) continuous mapping theorem (cf. Theorem 


3.4.4 in Whitt ( 20021 ). Rillingslev ( 19991 )). we have 


f^/nV^loo) 


9n 


h{u)du 


I 


V{oo) 


h{u)du. 


It is straightforward, from the change of variable, to check the left-hand side above af¬ 
ter taking an expectation equals to y/nP^ = ^/n]F[F^{V^{oo))]. From Theorem [T](b) 
and (A5’), we conclude that h{u)du : re > 1} is a uniformly integrable family 

and hence the first assertion follows. The second assertion in m is immediate from an 
asymptotic relationship between the scaled queue-l ength and offered waiti ng time processes, 
Q^jy/n — AF” => 0 as re —)• oo (see Theorem 6.1 in iReed and Ward! (j2008l )). combined with 
the weak convergence results in Theorem [T](a) and the converging together theorem. □ 
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4 Uniform Moment Estimates 


The main step in proving Theorem [T] is the tightness result of the family of steady-state 
distributions (indexed by n > 1) of the offered waiting time processes. The key to the 
desired tightness is to obtain uniform (in n) bounds for the moments of the steady-state 
distributions. Proposition [1] below provides estimates on moments of the (scaled) state 
process that are uniform in the scaling parameter n. 

In this section, it is helpful to emphasize the initial condition of the scaled process 
yn_ Pqj, initial condition Xq of in ([1]), y/rix^ —>■ a: as n —>■ oo, where x is the 
initial condition of the diffusion P in dSj). This is true if, for example, Xq = xl^/n for a 
given X > 0. Without loss of generality, and for notational convenience, we will assume 
Xq = xjy/n so that x is the initial condition of P”, and we will write P^. 

We follow the same convention for all processes defined in this section, and use the 
subscript x to remind the reader that the process definition assumes the initial state of Vff 
is X. 


Proposition 1. There exist Nq G IN and tQ G (0,oo) such that for all t > to 


lim sup — IE 


{v^{tx)y 


= 0 , 


( 12 ) 


where q > 2 is as in (A5), and q = 2 when (A5) is not assumed. 


We first prove Proposition [Hand Theorem [T] under the assumption that the hazard rate 
is scaled (that is, (A1)-(A4) hold for part (a) and (A5) additionally holds for part (b)). 
The case that the hazard rate is not scaled is covered by minor adjustments to the proofs, 
which are shown at the end of proof of Theorem [TJ 

The proof of Proposition [1] relies on a martingale representation of the offered waiting 
time process and the continuity properties of a nonlinear generalized regulator mapping. 
We first provide this setup, and second give the proof. 


Martingale Representation of the Offered Waiting Time Process 

Recall that we consider a sequence of systems indexed by n > 1 in which the arrival rates 
become large and service times small, and also the system primitives (t”, : j > l)n>i 

with t'j,v'j,d^ being the j-th arrival epoch, service time, and aban d onme nt time in the 


j ■ j ' j ^ “ * ' --- - 1 

n-th system. Our formulation closely follows that of iReed and Ward! (j2008l ). 


Dehne the u-fields where 


= u((^^ <),..., (tf, <, C+i) ^ 
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and let Tq = (T(t”)- Notice that is J^j?l_^-measurable and the abandonment time 

df of the i-th customer is independent of Hence, 


z = i, 2 ,..., 


(13) 


holds almost surely, where F” is the distribution function of d”. We then have a martingale 
with respect to the hltration (F”)i>i given by 


i=i 


Using (fT3l) . we also see that for alH G iV 


(14) 


Ml 


U(') = E 

i=i 


Next, consider the following centered quantities: 


«”(') ^ - !)• ^ ^ E(“-1)1 


{V"(t;;-)>dp- 


i=i i=i 

Algebra from ([I])-® and (fT3)) - (fT6]l shows that 

V:(t) = Xl(t) + e^(t) - ^ h-{u)du^ ds + Il{t), 

where 


(15) 


(16) 


(17) 


xlit) = ^+4 +(4) + (4 -1 

^/n n \ld / \ n / 


nt / rV^is-) \ 

e"(t) = J \J h'^{u)du\ds — 

II (t) = / l[yj*(s)=o]f)^'S. 

Jo 


1 

iiPI I 

Id Jo 


F^{V^{s-))dA-{sl 


Regulator Mapping Representation 

We define fluid-scaled and diffusion-scaled quantities to carry out our analysis. Let 

A^{t) = A^{t) = VIi(- -A’"A , and also 

n \n n 
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S^t) ^ M 2 ,At) = ^M 2 ,A[ntA 

Y n 

C(i) ^ v^y"(t), X2{t)^v^x2{t), ^At)^V^e2{t), I2it)^v^i2{t). 

The diffusion-scaled offered waiting time process can be written in terms of the fluid- 
and diffusion-scaled quantities as: 


v,At) = x2{t)+rAt)- 




t ( 


h{v)dv ds -|- 12{t), 


where 

X2{t) = x + 

^At) = r 


n 




s-{A^{t)) - s2,AAAt)) - M2,AAAt)) + 


+ 


n 


d" 


n 


'•VT(s-) 


h{v)dv ds — 


— [ —F^(—V"-(‘ 


dA^{s) 


It is straightforward to check using the definition of the idle time process that the process 
(V^'Ax) satisfies the conditions for the one -sided nonlinear generalized regulator mapping 
(cf. Definition 4.1 of Reed and WardI ( 2008I T {4>^, : D([0, oo),lR) —)• D([0, oo), [0, oo) x 

[0, oo)) is defined by {(jA = {z,l) where 


(^aj z{t) = y{t) — Jp ^ h(u)du'j ds + l{t) G [0, oo) for all t > 0; 
(b) I is nondecreasing, l{0) = 0, and z{t)dl{t) = 0. 

Hence 

{vaJ 2) = {c^\A’^){x2+A2). 


(18) 


n 


The proof of Proposition [T] is inspired by Budhiraia and Lee ( 20091 ). We expand their 
approach to incorporate customer abandonments. The key differences are (a) the use of 
the nonlinear generalized regulator mapping to connect results on the underlying arrival 
and service renewal processes to results on the offered waiting time process, (b) the use of 
basic properties of the underlying hazard rate function, and (c) a further consideration of 
the scaled error process e^. 


Proof of Proposition [T]. Let 

N2{t) = (l"(t) + 5"(^"(t))-5J,JA"(t))-M,"J^^ 




12 















so that X'^ can be succinctly written as 

x^{t) = x+b-t+N:^{t). 


(19) 


Recall that 6” converges to 6 */A as n —)■ oo from (A2). For the martingale in (fT^ . 
we note from Section 5.1 of Reed and Ward ( 20081 ) that both S2x{') and M2^{-) weakly 
converge to zero process as n —>■ oo, and also for all T > 0, lim^^oo sup^gjQ j’j |^”(^) — 
At| = 0 almost surely. These facts, together with the functional central limit theorem 
o An) ^ vai{ui)Bi, ■\/Avar(ni)R 2 ) imply that X" weakly converges to a 

Brownian motion with infinitesimal mean 6/X and variance cr^ = A var(Mi) + A“^ var(ni). 

For n > 1, fix a; G [0,oo) and write 14" (recall (fTSi) and (fT9l) l as 


R4(i) = ^^{x + b-i + N2{-) +?2(-) it), t > 0, 


( 20 ) 


where i : [0,oo) —[0,oo) is an identity map, i.e., i(t) = t for all t G [0,oo). Define a 
deterministic dynamical system analogous to (I20p as 


Z"(t) = (j)^ {x + h^i) (t), t > 0. 


( 21 ) 


Then, using the Lipschitz property of 4>^i-) (cf. Proposition 4.1 of Reed and Ward ( 20081 )1. 
we have for some k G (0, oo) that 


\Kit) - Zxit)\ < K sup |iV^(s) + e2(s)|, for all t > 0. 

0<s<t 


( 22 ) 


Next, denote by 


6 * ( 2 ) = 6" — / h{u)du for z > 0. 

Jo 


We claim that there exists a constant D > 0 such that Z^it) = 0 for all t > Dx. Notice 
that ^(&n('Z)) = —hiz) < 0, then from (A4) we conclude that there exist z, 5 G (0, 00 ) and 
Ni £ IN satisfying 

inf inf |6;(z)l > (5. (23) 

The above condition ( 1231) resembles the one-dimensional version of the so-called ‘cone 
condition’ 


m 


Atar et al. (200 il l that was used to characterize stability of a family of 


diffusion models with state dependent coefficients, constrained to take values in some 
convex polyhedral cone in The condition (I23p says that the state-dependent drift 

(6 *(z),n > Ni,z > z) stays away from zero (to the negative direction), uniformly both in 
the parameter n and the state variable z. Thus, for all n > Ni and z > z, 


6* (z) G €{ 5 ) = {n G [0, 00 ) : |n| > <5}. 


(24) 
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For go £ [0, oo), denote by T(go) the collection of all trajectories of the form 
4>{t) = Qo + [ C{s)ds + l{t) £ [0,oo), for all t > 0, 

Jo 

where l{-) is nondecreasing, 1(0) = 0, 4’{t)dl{t) = 0, and ^ : [0, oo) —)■ iR is a measurable 

map satisfying for all t € [0, oo), |^(s)|ds < oo, and ^(t) € C(S). Define the hitting time 
to the origin function 

T(go) = sup inf{t > 0 : (/>(t) = 0}. 

<^eT(go) 


Then, Lemma 3.1 of lAtar et ahl (j200lh (in particular, the estimate (3.10) therein) shows 
that 

T{qo) < and for all (p G T{qo), 4>{t) = 0 for all t > T(go), 


where k G (0,oo) is as in (12211 . Combining this observation with (I24p . we now see that 
Z^{t) = 0 for all n > A^i, t > Dx, where D = jd. Using this estimate in ([22]) . we now 
have that 

\V^(tx)\<K sup |iV,"(.)+?2(.)| (25) 

0<5<tir 

for all re > A^i, t > D and for all initial conditions x G [0,oo). 

Next, we obtain an estimate on the g-th moment (g > 2) of the right-hand side of (1251) . 
For the martingale in (|25p . it follows from the discussion below (I19p that there exists 
N £ IN such that for all n> N, 


IE 


sup |iV;(s)|^ 

0<s<t 


< co(l + t) for some cq £ (0, oo). 


(26) 


(Note that such an estimate is available under the second moment condition lE{ui + vi)‘^ < 
oo.) Moreover, under the g-th moment (g > 2) assumption on the primitives as in 
(A5). we have from sta.ndard estimates for renewal processes (see, e.g.. Theorem 4 in 
Krichagina and Taksar ( 1992l U and observations leading to (1261) that there exists N £ JN 
such that for all n > N, 


IE 


sup |iV;(s)r 
.0<s<t 


< ci(l + for some ci G (0, oo). 


(27) 


Also, since ^(O = 4> 0 uniformly on compact sets as re —)• oo (see Section 5.1.2 of 
Reed and Ward ( 2008l l. and note that our assumed non-zero initial condition is immaterial), 
there exists an N 2 £ IN such that for all re > N 2 , 


IE 


sup i?s(s)r 

0<s<t 


< C 2 for some C 2 £ (0,00). 


(28) 


14 

























Applying the estimates in (l26]l . (1771) and (l28l) . we now have that for al\ q > 2, t > D and 
X G [0,oo), and n > Nq = max{N, Ni, N 2 }, 

lE\V^{tx)\'^ < 03(1 + for some C3 G (0,00). (29) 

Then the desired result ()12p follows on choosing to = D. □ 


5 The Steady-State Convergence Proof 


We begin by providing a general statement concerning Markov processes which we require 
to apply in the steady-state convergence proof. For S G (0,oo), define the return time to a 
compact set C C [0,00) by = inf{t > 6 : V^{t) G C}. 


Dai and Mevn 


Prop osition 2. fTheorem 3.5 of Budhiraia and Lee l200d) . cf. Proposition 5-4 
i '199,i) ) Let f : [ 0 ,oo) —>■ [ 0 ,00) be a measurable map. Define for 5 G ( 0 ,00), and a compact 
set C C [ 0 ,00) 

r rcCs) .. 

Gn{x) = lE / f{Vfi{t))dt , xG[0,oo). 

Jo 


If sup„ Gn is everywhere finite and uniformly bounded on G, then there exists a constant 
rj G ( 0 ,00) such that for all n G IN, t G [ 5 ,00), x G [ 0 ,00), 


-]E[G, 


.(r"(())l + i f ]Elf{v"{s))]d 
^ Jo 


s < -Gn{x) + T]. 


(30) 


Proof of Theorem [H Part (a): Let tt”’ and vr be the unique steady-state distributions 
of P” and P in dHI) respectiv ely. Since the process- level convergence P"" => P as n —)• 00 
follows from Theorem 5.1 in Reed and Ward ( 20081 ). it suffices to establish the tightness 
of the family {vr"'}. To prove that tightness, it suffices to show that there exists a positive 
integer N such that for all re > V 



X7r'^{dx) < c, 


( 31 ) 


where c G ( 0 ,00) is a constant independent of re. 

A key observation from Proposition [ 1 ] is that there exists a 70 G ( 0 ,00) such that, for 
to and No as in that same proposition. 


sup lE{Vfi{tox)y < \x^^ for x G (70,00) and q = 2. ( 32 ) 

n>No 2 
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Next, we intend to apply Proposition [2] with 

5 = to 7 o, f{x) = 1 + for x £ [0, oo), and C = {x G [0, oo) : x < 70 }, (33) 

where q = 2. (In the proof of part (b), q will be greater than 2.) 

Suppose we can show that there exists a positive integer N and c G (0, 00 ) such that 


sup E 

n>N 


/ (1 + (C(t))'?-i)dt 

Jo 


<c(l + x'^), xG[0,oo), 


(34) 


so that the conditions of Proposition [2] are satisfied. Then, for x G [0,oo) and t] G (0,oo) 
the constant in that same proposition, 

chjx) = \Gn{x) - \E[Gn{Km >\f^ E{f{V^{s)))ds - V, 

for all n G JN, t G [d, 00), and x G [ 0 ,00), and so 


/ 4>n(x)7r"-((ix) > / 

'[OjCxd) J[0,oo) 


jJ^Eif{V,-is)))ds-ri^7r-{dx). (35) 


Furthermore, it follows from the definitions of a steady-state distribution vr" and a function 
<h„(x) that 0 = JjQ $„(x) 7 r”(dx) if is bounded and measurable, and, otherwise, 


Fatou’s lemma shows that 


0 > / $„(x)7r”(dx). 

J [0,oo) 


(36) 


Fubini’s theorem and the definition of a steady-state distribution show 


^ f E{f{V^{s)))dsi:-{dx) = 7 r / E{f{V^{s)))ir-{dx)ds = [ /(x)vr"(dx). 

70 , 00 ) ^ Jo ^ Jo ./[ 0 ,oo) ./[Ojoo) 

(37) 

Finally, it follows from (I5^ - ([571) that 


0 > 


70 , 00 ) 


/(x) 7 r”(dx) -7, 


(38) 


which establishes (1311) . 

It remains to show ()34[) . Gi ven the estimate in (1321). t he proof follows along the same 


lines as those of Theorem 3.4 in Budhiraia and Lee ( 2009li . because <5, f{x) and G in ([331) 


are chosen in the sam e way as in the cited the orem. Following steps analogous to those 
leading to (39)-(41) in iBudhiraia and Led (120091 ). it suffices to check that the deterministic 
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trajectory in m satisfies the following property: For every ci G (0,oo), there exist 
N £ IN and C 2 G (0, oo) such that 

sup sup |Z”(t)| < C 2 (l + x), xG[0,oo). (39) 

n>N 0<t<cix 

From the Lipschitz continuity of the generalized regulator map cj)^, we have 


sup \Z^{t)\ < k{x + sup for some fi:G(0,oo). 

0<t<cix 0<t<cix 


Letting C 3 = sup„>;^ | 6 ”| for sufficiently large N G FV, we have \b'^t\ < c^t since ^ 6/X £ 
(— 00 , 00 ) as n —>• 00 . Then, the desired estimate (l3^ follows on setting cy = k + C 1 C 3 . 
Therefore, the key step (41) in the proof of Theorem 3.4 Budhiraia and Lee ( 2009l f is now 
satisfied and this completes the proof. 


Part (b): Under (A5), Proposition [T] implies (|32l) with q > 2. Therefore, we have a uniform 
(in n) moment bound on 1E[V ^obtained in (|38p . with g > 2 as in (A5). Combining 
the weak convergence result in part (a) and a uniform integrability of {V^(oo)'^“^}„>i, we 
conclude the desired moment convergence result. 


The proof modifications when the hazard rate is not scaled are: 


• In the primitives used in Section 01 to obtain the martingale representation of the 
offered waiting time process, and are replaced by dj and F. 

• The dehnition of e” is replaced by 


• The diffusion-scaled offered waiting time process is now represented as: 

v:{t) = x2{t) +Fi{t) - fF'm:{s-)ds+ 

• The function h{u) is replaced by the constant A'(0). Then the nonlinear regulator 
mapping cj)^ is going to become a linear regulator mapping, and the definition of 6 * 
is going to change so that the second term is linear. 


□ 
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